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Abstract 

In [5] we studied an interacting particle system which can be also 
interpreted as a stochastic growth model. This model belongs to the 
anisotropic KPZ class in 2 + 1 dimensions. In this paper we present 
the results that are relevant from the perspective of stochastic growth 
models, in particular: (a) the surface fluctuations are asymptotically 
Gaussian on a Vint scale and (b) the correlation structure of the 
surface is asymptotically given by the massless field. 

1 Introduction 

We consider growth models of an interface described, at least on a coarse- 
grained level, by a height function x i— >■ h{x,t), where t G M is the time and 
X G M'^ (or Z*^) is the position on a d-dimensional substrate. When the growth 
is local, then on a macroscopic level the speed of growth f is a function of 
the slope only 

v = v{\/h). (1) 

The class of models we consider are stochastic with a smoothing mechanism 
leading to a non-random limit shape, 

hUO ■■= lim (2) 
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and in the Kardar-Parisi-Zhang (KPZ) class, related to the KPZ equation 



^ = uAh + ^iVhr + v, (3) 

with 7] a local space-time noise. ([3]) is a mesoscopic equation, in which the 
smoothing mechanism is modeled by the uAh term, which is physically re- 
lated with the surface tension, t] is the stochastic part, while the non-linear 
term |(V/?.)^ is coming from Taylor series of v(Vh) around V/i = 0. Indeed, 
([3]) assumes that one sets the system of coordinates such that V/i is small. It 
is then possible to make the choice so that the 0th and 1st Taylor coefficients 
vanish. In [26] it was argued that for the long time behavior, higher order 
terms become irrelevant. 

Before describing the two-dimensional case, let us spend one paragraph 
on the d = 1 situation. By non-rigorous arguments, the scaling exponent 
in d = 1 were correctly derived to be 1/3 for the fluctuations and 2/3 for 
the spatial correlations, see [3,4,20]. This means that the height function h 
rescaled as 

, resc. ^ _ /^(^t + Ut^-^ t) - th^, ((et + Ut'/-')/t) 

in the long time limit will converge to a well-defined stochastic process. The 
analysis of simplified models in the KPZ class gives the following results for 
noise-free initial conditions (results expected to hold for the whole class by 
universality). If /ima is straight (i.e. second derivative is zero) around ^, then 
the limit process is the Airyi process [7, 8] , if /ima is curved, the process is 
the Airy2 process [23,24,36,38]. Further results as the transition between 
these processes [10] or correlations at different times [6,9,22] with also some 
unexpected results [15] are also available. See [12, 16, 17, 34] for reviews 
and [1,2,18,25,32,35,37] for related works. 

The main reason of this paper is to present results on the anisotropic 
class in d = 2. Going back to the KPZ equation, one realizes that in two 
dimensions one should write 

dh 

— = i^Ah + Q{Vh)+r], (5) 

where Q is a quadratic form. It turns out that there are two classes, depend- 
ing on the signature of the quadratic form Q. 

(a) Isotropic KPZ: sign(Q) = (+, +) or (—,—), i.e., the Hessian of 
V = v{Wh) has two eigenvalues of the same sign. In this case, the fluctuations 
grow as for some /5 > 0, but renormalization arguments used in the d = 1 



2 



Figure 1: A configuration of tlie model analyzed with = 100 particles at 
time t = 25. In [14] there is a Java animation of the model. 



case fail to give an answer. Numerical studies on discrete models belonging to 
the KPZ class [19,40] as well as numerical solution of the equation itself [31], 
indicate that (3 = 0.240 ± 0.001 (ruling out the value 1/4). For this case, 
essentially no analytic results are available. 

(b) Anisotropic KPZ (=AKPZ): sign(Q) = (+, — ). In this situation, 
the behavior is quite different and the fluctuations increase much slower, 
only as Vint. The first time the AKPZ was considered was by Wolf in [41], 
who was interested in growth on vicinal surfaces. They are surfaces with a 
small tilt with respect to a high-symmetry plane of the crystalline structure 
of the solid, and they are intrinsically asymmetric. By detailed one-loop 
computations. Wolf deduced that the non-linearity should be irrelevant for 
the roughness, which should grow in time as ln(t), i.e., 

Var(/i(^t, t)) ~ ln(t), as t cx), (6) 

exactly as in the A = case (Edwards- Wilkinson model [13]). On the nu- 
merical side, this drastic difference have been tested positively soon after [21] 
on the model used by Wolf directly, but also on a deposition model [30]. On 
the analytic side, Prahofer and Spohn [33] considered a microscopic model 
in the AKPZ class and reproduced Wolf's prediction, namely ([6]). Moreover, 
they computed the local correlations obtaining 

lim Var(/i(x, t) — h{x' , t)) ~ In |x — (7) 

t— >oo 

for large \x — x'\ —>■ oo, but not growing with t. 
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Figure 2: Initial condition. The corner of the crystal is at position (0,0,0). 



In our recent work [5] we consider a growth model in the AKPZ class, 
for which we determine the detailed correlation structure on a macroscopic 
scale. Below we explain the model and present the results relevant under the 
perspective of growth models (see Figure [T] for an illustration of the surface 
we analyzed). 

2 The model 

We consider a simplified model of anisotropic evaporation of a crystal. As 
initial condition we have a perfect corner on the positive octant, see FigureO 
The dynamics is as follows. The evaporation occurs by detaching of columns 
along the X2-direction. Each cube at {xi,X2,X3) such that both its X3-facet 
and Xi-facet are visible, evaporates with rate one, and by doing so it takes 
away also all the cubes in the X2-column above it: {{xi , X2,X3) S.t. X2 ^ '^2}' 
This dynamics is intrinsically asymmetric, since the X2 and X3 directions play 
different roles. 

There are several possible choices of coordinate which can be used to 
describe the interface as in Figure [TJ The one we have chosen in [5] is the 
following. Consider the two-dimensional projection as in Figure [H whose 
points can be described by (xi — y,y = xs + X2). For convenience, we take 
X = xi as the first coordinate and y as the second one, while the height 
function will be given by 0:3: 

{x, y, ht{x, y)) = {xi,X2 + X3,X3). (8) 

The surface grows with a non-zero speed. Thus we rescale our coordinates as 
Xi = C,t, X2 = Cty ^3 = ht{C,t, (t) so that X = ^t, y = (t — ht{C,t, (t). Moreover, 



4 




Figure 3: An example of the first three moves. The moved cubes are the ones 
with the black circles, and also the ones with white circles which are moved 
because they are in a column above the black ones. 



in the long time limit, our surface has a deterministic limit shape 

ftK,C):=li^MI^. (9) 

The hmit shape consists of facets joined by a curved region in the coordinate 
range 

{i-vcr<^<{i+vcf- (10) 

Explicitly, for (^, Q in the curved region, 

Me, C) = - arccos 2v/C J ^^^^ 

In the curved region, in the long time limit, one locally sees a flat interface 
characterized by the two slopes 

dh dh 

The two slopes determine imiquely a stationary measure, and since our growth 
is short range, the speed of growth v will be a function of and only. 
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Indeed, we get 



1 sin^Tih^) sin(7r/i^) 



(13) 



From this expression, we can compute the determinant of the Hessian of v 
with the result 



-An' 



< 



(14) 



for h^, h(^, h^ + h^ G (0, 1) (which is equivalent to being in the curved region). 

Thus, our stochastic growth model has a smoothing mechanism determin- 
ing deterministic limit shape, the growth rule is local, and the Hessian has 
eigenvalues of opposite signs. Therefore, our model belongs to the anisotropic 
KPZ class (AKPZ) in 2 + 1 dimensions. 



3 Results 

Logarithmic fluctuations 

The first result concerns Wolf's prediction of the logarithmic growth of the 
fluctuations. In Theorem 1.2 of [5] we prove the following result. 

Theorem 1. Consider any macroscopic point in the curved region, i.e., 
(x, y) = {C,t, (t) with (^, C) satisfying flU) . Then, the random variable 

r (c /-^ ht{x,y) - E(/j-f(x,y)) 

ht{iX)-= 7f= (15) 

Vint 

converges in the t —>■ oo limit to a normal random variable with mean zero 
and variance l/(27r^). 



Complex structure and local measure 

To present the other results, we need to define the following mapping. For 
any () satisfying ffTUl) . define the map Q to the upper (complex) half-plane 
M.:= {z eC s.t. lm{z) > 0} as follows: 

\n\ = y^, \l-n\ = y^. (16) 

Geometrically, the point is the intersection in H of the circle^] B{0, ^/C) 
and B{l,^/^) (see Figure Hj). Also, denote the angles of the triangle with 

"'^The notation B{x, r) denotes the circle centered at x of radius r. 
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vertices 0, 1, by 

TT^ = — arg(l — Q), vr^ = arg(fi), tti = vr — tt^ — tt^. (17) 

Now we discuss the local measure in the curved part of the surface. If one 
focuses around macroscopic point (), asymptotically one sees a translation 
invariant measure parameterized by the gradient of the surface. The surface 
can be also viewed as a random tiling of the plane by three types of lozenges 
indicated in Figure HI Then, the measure one sees locally is the unique trans- 
lation invariant Gibbs measure on lozenge tilings with prescribed proportions 
of lozenges (vr^/vr, vr^/vr, tti/tt) (see [11,27,29] and references therein for dis- 
cussions of these Gibbs measures). From the tiling point of view, it is easy 
to understand the connection between the angles tt^, tt^ and the asymptotic 
slopes h^, h(^: 

= — vr^/vr, /i^ = 1 — vr^/vr. (18) 

Indeed, along the ^-direction, the height changes whenever there are lozenges 
associated to vr^, while along the ( direction it changes whenever there are 
lozenges not associated with vr^ (varying ( with fixed ^ corresponds to moving 
along the (1,-1) direction in our figures). Finally, the growth velocity takes 
a very simple form in terms of Q: 

v = lm{Q)/n. (19) 

Interestingly enough, although we study a stochastic growth model (out 
of equilibrium), the limit shape is related to a variational problem appearing 
for uniformly distributed lozenge tilings of a given domain. The map Q 
satisfies the complex Burgers equation 



which has been identified as the Euler-Lagrange equation for the surface free 
energy for fixed domains, see [28] and references therein. The hmit shape 
gives a solution of ( l20l) and it is encoded in the equation 

^n + cii-^) = ^ii-^), i^gh, ^,c>o. (21) 

Indeed, for any i7 G H, there is a unique couple (^, C) ^ ^+ satisfying ( l2Ti) . 
and then the limit shape h{^, () can be simply expressed as the integral of 
the slope as in (|TT1) . 

Macroscopic correlations 

Finally, let us discuss the covariance structure of different macroscopic points 
in the curved part of the surface. At first sight, in view of ffT^ one might 
think that we have to compute E(nj'!^-^ ht{C,i, Q)) for distinct Q) satisfying 
ffTUj) . However, this turns out to converge to zero in the t —>■ oo limit. 

Our result is that, after mapping by Q, the height function converges to 
the massless field on H, which we denote GFF. The value of GFF at a point 
cannot be defined. However, one can think of expectations of products of 
values of GFF at different points as being finite and equal to 

E[GFF(^i)---GFF(^J] 

{0 if m is odd, 

pairings a 

where 

E[GFF(2) GFF(w)] = 

is the Green function of the Laplace operator on EI with Dirichlet boundary 
conditions. In fl22|) one recognizes Wick formula for Gaussian processes. For 
further details on the massless field (also called Gaussian Free Field), see 
e.g. [39]. 

Our result on the correlations on a macroscopic scale is as follows (The- 
orem 1.3 in [5]). 

Theorem 2. For any m = 1,2, . . ., consider any m distinct pairs Ci),i = 
1, . . . , m} which satisfy / flOj) . Denote 

Ht{^, C) := M^t, Ct) - Eihti^t, Ct))] (24) 

and Qi := Q{C,i,(i). Then, 

lim E{Ht{^i, Ci) ■ ■ ■ HtiU Cm)) = nGFF{n,) ■ ■ ■ GFF(a^)] (25) 



Giz^{m-i),z^{m)) if m is even, (22) 



Q{z,w) :-- 



27T 



z — w 



z — w 



(23) 
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At first sight it might look strange that in Theorem [T] we have vhit 
scahng, while in Theorem [2] it is not present. If we consider Ht as random 
variable, Theorem [T] tells us that, in the t —* oo limit, Ht does not converge 
to a well-defined random variable (the Vln/f: normalization is missing). On 
the other hand. Theorem [2] tells us that the macroscopic correlations of Ht 
are non-trivial. The reason is that asymptotically the interface Ht will take 
values in the space of distributions, not ordinary functions. 

Our Theorem 1.3 in [5] is actually more general. Indeed, we obtain the 
limiting correlations also for the height function at different times, provided 
that the space-time events at which we consider the height function satisfy 
a "space-like condition". When we consider heights at different times, we 
cannot use t anymore as a large parameter. We use L instead, set t = tL, 
and denote by [vL^riL^TL) a space-time point (z/ = r,^, r] = tQ. Then, the 
condition ( fTOl) that a point is in the curved region becomes 

iV^-Vv)' <'^<{V^+Vv)'- (26) 

Then, the mapping Q is now from a triple (z/, rj, r) to EI given by the condition 

\n\ = ,/^, \i-n\ = y^. (27) 

Then the full result of our Theorem 1.3 in [5] is the following. 

Theorem 3. For any m = 1,2,..., consider any m distinct triples 
{{ui, rii,ri),i = 1, . . . , m}, which satisfy / f^) and the space-like condition 

Tl < T2 < . . . < T„„ r]i>r]2>--->Vm- (28) 

Denote 

Hl{u, t], r) := [Kl{i^L, r]L) - E(/i,i(z^L, r]L))] (29) 
and Qi := Vtiyi^rji^Ti). Then, 

lim E(i7i(z/i, r/i, n) ■ ■ ■ i^i(^^„^, Vm, r^)) = E[GFF(^]l) ■ ■ ■ GFF(fi„)]. (30) 

L— >oo 

There is an interesting fact to notice. Consider for simplicity the following 
two pairs of space-time points 

{(i^i, Vu 1), (i^2, V2, 1) with ri2 < rji} (31) 

and 

{(z/i,r/i,l),(z/',r7',r) s.t. ^{1^2,^2^) = ^(i^' ,v' ,r),T > 1}. (32) 

Then, Theorem [3] tell us that the correlations of Hl at these two pairs are 
asymptotically the same. Also, remark that space-time points have the same 
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image by 1) on H if and only if they have the same macroscopic slopes defined 
in dn]), see (dHD- 

This a priori unexpected phenomenon, can be explained as follows. Along 
the space-time curves with constant Q, the height decorrelates on a much 
slower time scale than L. The same fact has been already observed in KPZ 
growth on a (i = 1 substrate [15] (there the correct scale is L instead of L^/^ 
for spatial correlations). 

As a consequence of this interpretation, we conjecture that Theorem [3] 
holds without the space-like condition ( |28l) . provided that Qi, . . . , Qm remain 
pairwise distinct. 
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